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Abstract: Jet emission at hadron colliders follows simple scaling patterns. Based on per- 
turbative QCD we derive Poisson and staircase scaling for final state as well as initial state 
radiation. Parton density effects enhance staircase scaling at low multiplicities. We pro- 
pose experimental tests of our theoretical findings in Z+jets and QCD gap jets production, 
based on minor additions to current LHC analyses. 
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1 Multi-jet rates 



Multi-jets final states are ubiquitous at hadron colliders. QCD jet radiation mostly off 
the initial state partons has huge impact on almost every LHC analysis. For example 
in Higgs searches an accurate description of the jet recoil against the Higgs boson allows 
for an efficient rejection of many backgrounds [1]. In top pair production and single top 
production the identification of QCD recoil jets on the one hand and decay jets on the other 
hand is one of the limiting factors in precision analyses [2, 3]. New physics searches largely 
rely on hard decay jets of new strongly interacting particles, which makes them vulnerable 
in the case of soft decay jets hidden in the QCD jet activity [4-6]. Understanding the jet 
multiplicity and the jet spectra from QCD radiation is a core ingredient to improving any 
of these analyses. 

Even though jet radiation seems to follow simple patterns [7, 8], theoretical predictions 
for multi-jet observables in perturbative QCD are challenging. It is clear that the radiation 
of relatively hard quarks and gluons is a direct consequence of the hierarchy between 
the large proton-proton collider energy and the typical electroweak scale of the partonic 
interaction [9, 10]. Numerically, we can combine the QCD parton shower with hard matrix 
element calculations, to predict jet radiation patterns over wide phase space regions [11— 
13]. Analytically, Sudakov factors and generating functionals can be used to describe 
QCD jet radiation [14]. A careful comparison of analytical and numerical approaches to 
LHC data would allow us to determine the strengths and limitations of the underlying 
theoretical concepts: fixed-order perturbation theory, parton showers, and resummation 
based on generating functionals. 

In fixed order perturbation theory leading-order jet rates are available for effectively 
arbitrarily high multiplicities [15, 16]. The number of jets we can consider is limited only by 
computing power. However, leading order predictions suffer significant shortcomings when 
it comes to precision. The renormalization scale dependence as a measure for the theoretical 
uncertainty grows with each power of a s when we add a final state jet. High powers of the 
logarithmic scale dependence (log (j,r) u mean that such cross section predictions can only 
be considered an order-of-magnitude estimate. In addition, in the presence of phase space 
constraints large logarithms spoil the convergence of fixed order perturbation theory. For 
example jet vetoes will induce sizeable logarithms of ratios of leading jet px or masses (mj? 
or mz) to additional resolved jets (experimentally relevant down to 20 — 30 GeV) [17, 18]. 

Next-to-leading order computations ameliorate the scale dependence and capture an 
additional logarithm. However, we are still limited in the available final-state multiplicity, 
e.g. pure jets are available for nj e ts < 4 [19], nj ets < 4 in association with W/Z bosons [20, 
21], for ti production and Higgs production in gluon fusion NLO corrections are known for 
Wjets < 2 [22, 23]. However, over the past few years this field has progressed enormously. 
As a consequence for Standard Model processes a similar level of automation as for leading- 
order calculations is within reach [24]. An approach applicable for general New Physics 
extensions, though limited to 2 — > 2 processes, has been presented in [25]. 
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Although NLO calculations contain one additional power of enhanced logarithms, this 
might not be sufficient for high jet multiplicities. At NNLO, although there has been an 
enormous amount of recent development [26], the number of fully differential calculations 
is limited, and an automated implementation is not foreseeable in the near future. 

On the other hand, we know that jet radiation is enhanced by traceable logarithms. 
This makes improved predictions for QCD observables based on resummation possible. 
The general strategy is to redefine the perturbative series from powers of a s to including 
the relevant logarithms; the simplified structure of these enhanced terms then allows for 
a resummation to all orders. Once the resummed form is known we can match onto a 
fixed order calculation and avoid double-counting. For Sudakov-type logarithms a general 
method for this type of resummation is available [27], and for particular event shape ob- 
servables an automated approach exists [28]. In LHC analyses, the resummation of finite 
logarithms in the presence of a jet-veto scale is of interest [17, 29]. 

A numerical approach to resummation is provided by parton-shower simulations [30]. 
It is automated in the multi-purpose Monte Carlo generators Pythia [31], Herwig [32] and 
Sherpa [33] to leading order in the strong coupling combined with the resummation of lead- 
ing collinear logarithms (LO/LL). This method differs from the previous approaches in that 
the full spectrum of final state partons or hadrons is produced explicitly. While the parton 
shower is well defined for relatively small transverse momenta of the jets it is not applicable 
for hard jet radiation. However, this limitation is overcome by the CKKW [11], MLM [12], 
and CKKW-L [34] jet-merging algorithms, that incorporate the tree- level matrix-element 
corrections for the first few hardest emissions [30, 35]. 

A complementary strategy is provided by the MC@NLO [36] and POWHEG [37] ap- 
proaches, that realize the matching of NLO calculations with parton showers. While these 
methods guarantee NLO/LL accuracy only the first /hardest shower emission gets corrected 
by the real-emission matrix element. Higher jet multiplicities are described in the parton- 
shower approximation only. First attempts to combine the NLO/LL approaches with the 
tree-level merging ansatz have been reported recently [38] . An unprecedented level of so- 
phistication for predicting multi-jet final states is achieved by the promotion of merging 
algorithms to next-to- leading order accuracy [39]. 

Even though we can nowadays simulate multi-jet events, a detailed understanding of 
inclusive or exclusive nj e t s distributions at the LHC is still missing. Its universal features 
have been studied since 1985 [7]. Scaling patterns can be conveniently displayed in the 
ratio of successive exclusive jet cross-sections 

R(n+l)/n = = W1 t h Pn = • (1-1) 

We define the jet multiplicity n as the number of jets in addition to the hard process, 
e.g. o\ for pure QCD di-jets is experimentally a 3-jet final state. Jets which are part of the 
hard process are not included in the scaling analysis because they do not arise from single 
QCD emissions. 
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Two patterns provide limiting cases for most LHC processes and are referred to as 
staircase and Poisson scaling. Staircase scaling is defined as constant ratios between the 
successive multiplicity cross sections 

R(n+l)/n = R = e ~ b > (I- 2 ) 

where R and b are constant. The exclusive n-jet rate for this distribution is a n = aoe~ bn 
where do is the 0-jet exclusive cross section. Staircase scaling for exclusive and inclusive jet 
rates is equivalent, with identical values of R. For a Poisson distribution with expectation 
value n the rates are 

Pn = j— Or R(n+l)/n = ~r ■ (1-3) 

More properties of the distributions are described in Ref. [1]. 

Examples of Poisson [40] and staircase [41] scaling are found in current LHC analyses, 
although especially with the former, the examples are limited. Recent computations in 
fixed order perturbation theory [42] confirm that the staircase pattern improves at NLO 
compared to the LO prediction. In this paper we propose a number of different test-beds 
for both scaling behaviors, and provide predictions for the nj e t s ratios. In particular we 
suggest two explicit measurements based on simple extensions of already existing analyses, 
namely Z+jets production and gap jet studies in pure QCD events. 

Establishing the origin of scaling in radiated jets from first principles is the primary 
purpose of this publication. The line of reasoning was touched on in our previous work, but 
a complete treatment is lacking. Besides constituting an aspect of perturbative QCD on 
its own, a fundamental understanding of scaling is obligatory for subsequent phenomeno- 
logical applications. Only through an underlying mechanism can we convincingly argue 
that scaling patterns supplement the information provided by fixed-order or parton-shower 
calculations. 



This paper is arranged as follows. In Sec. 2 we introduce the standard theoretical 
framework for the parton shower, and elaborate on primary emissions with respect to the 
core process as a necessary condition for Poisson scaling. In conjunction, we demonstrate 
in the context of generating functionals that scaling patterns naturally emerge and can be 
proved to all multiplicities. In Sec. 3 we generalize our discussion to hadron colliders. We 
introduce the corresponding generating functional formalism and discuss the effect of initial 
state PDFs, both through the kinematics and initial-state backward evolution. We close 
this section by showing that also in BFKL evolution the same scaling patterns appear as 
limiting cases. Having a firm understanding of the emergence of the scaling patterns from 
QCD, we can define cut regimes which exhibit idealized scaling. In Sec. 4 we propose two 
experimental measurements where the underlying hypotheses of this paper can be tested 
at the LHC. We conclude and discuss possible extensions of our ideas in Sec. 5. 
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2 Final state splittings 

Poisson statistics describe individual events occurring repeatedly and independently. For 
particle emissions this implies a non-trivial assumption. For successive photon radiation 
in QED it is the standard solution to divergences arising from soft and collinear photon 
radiation. The underlying approximation for multiple soft photon radiation is the eikonal 
approximation, in which the radiation of a soft photon does not affect the hard process. 
The eikonal approximation states nothing but a basic assumption of the Poisson process — 
each successive photon sees the same unaltered hard process. For QCD it has been known 
for some time that the non-abelian nature of soft-gluon radiation leads to a deviation from 
Poisson statistics for gluon multiplicities [43]. When we measure jet ratios it is important to 
separate general exponentiation [44, 45] from a Poisson distribution; to generate a Poisson 
distribution the exponentiation of soft real and virtual corrections has to include exactly 
the single-emission matrix element. 

2.1 Final-state parton shower 

The simulation of a LO/LL event in the parton shower approximation starts with the 
generation of a single phase space point for the partonic core process. The process' external 
(colored) lines then act as seeds for the subsequent parton shower evolution. Driven by 
unitarity they start at the hard process scale Q and finish at the shower cutoff scale Qq ~ 
1 GeV. Hard matrix element corrections for an arbitrary number of additional particles can 
be added to the parton shower using the above-mentioned matching schemes [11, 12, 34]. 
However, in this section we treat all additional emissions as coming from the parton shower 
and disregard matrix element corrections. 

The basis of the LL parton shower is the fully factorized form of the collinear matrix 
element and phase space da n+ \ ~ da n P\-+% d<§>\. Using this simplification the parton 
shower remains local, but loses information on spin, color correlations and interference 
effects in addition to higher order terms neglected in the 1 — > 2 splitting kernels Pi_>2- 
Besides practicality, one of the benefits of collinear factorization is that the resummation 
of LL and some NLL contributions follow very naturally. To see this, we represent the 
evolution along an individual line by integrating the Sudakov factor over the appropriate 
virtuality scales from the lower cutoff Qq to a free hard scale t 



For the regularized splitting kernels r,- we use the next-to-leading logarithmic approxima- 
tion 



with color factors Cj = Cp {Ca) and the constant terms Aj = 3/2 (11/6) for gluon emission 
off a quark (gluon). The lower cutoff scale Qq is omitted in the argument of the Sudakov 
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Figure 1. Simplest primary (left) and secondary contributions (right) assuming a core process 
with a hard quark line. 

form-factor. Expanding the exponential we see that Eq. (2.1) represents an arbitrary 
number of soft and collinearly enhanced emissions, either resolved or unresolved. 

To describe a parton-shower simulated event we note that the QCD evolution proceeds 
as an integration of the product Sudakov along the virtuality t, 

A(t)= [] A 3 (t)^e~ T . (2.3) 

ext lines 

The product defining T is over the appropriate factors for each external line, where j 
denotes the particle flavor. Limiting ourselves to final state splittings this expression only 
contains evolution kernels as shown in Eq. (2.1), and it is by construction guaranteed to 
exponentiate with an appropriate expression T. As long as T is fully local and does not 
depend on previous emissions it is guaranteed to produce a Poisson distribution for the 
multiplicities. The exponentiated form in Eq. (2.3) immediately identifies n = T. This 
statement does not depend on the form of T or its dependence on the hard scale t. All that 
matters is that each splitting does not change the subsequent evolution. In the remainder 
of this paper we define all emissions directly contained in the expansion of Eq. (2.3) as 
primary with respect to the core process. 

The first splitting in the parton shower picture defines the single emission probability. 
Following Fig. 1 a second emission can then appear from the original leg or off the first 
emission. For the former, this emission is contained in Eq. (2.3) and does not change 
the Poisson pattern. The latter changes the exponential; we refer to it as secondary with 
respect to the original hard process. From a scaling perspective the relevant questions are 
first, what is the relative size of the two contributions; and second if we can change the 
individual strengths of primary and secondary emissions through kinematic cuts. 

In the parton shower approximation we can associate specific integrals over virtuality 
with individual partonic structures appearing in the final state evolution. An alternative 
evolution ordered in a consistent variable (e.g. angle) is logarithmically equivalent. Using 
this formalism the primary contribution to two gluon emission off a hard quark shown in 
Fig. 1 is 

rQ 2 rQ 2 
(7 pr 1 maxy ( g2 ) g2 } = primary / ^ T(Q 2 ,t)A g (t) / d£ T{Q 2 ,t')A g (t') . (2.4) 

JQl JQl 

The coefficient c primary which includes the Sudakovs associated with the hard line is process 
dependent, as this hard line can be either a quark or a gluon. The two external scales are 
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the scale Q of the hard process and the lower cutoff scale Qo- If the primary emissions 
are strongly ordered in the evolution variable, the corresponding phase space factor 1/2 is 
absorbed in c P nmar y. The simplest secondary contribution also shown in Fig. 1 is, 

rQ 2 rt 
a secondary (Q 2 ,Qg) = c secondary / dtT(Q 2 ,t)A g (t) dtT(t,t')A g (t') . (2.5) 

JQl JQl 

The splitting kernels in the two expressions only differ in the integral boundaries for the 
second emission. In the leading logarithmic approximation (in the exponent) for the Su- 
dakov factors, we can perform the integrals in Eq. (2.4) and Eq. (2.5) in terms of error 
functions. The full expressions are not particularly enlightening, but two specific limits 
contain crucial information. 

(1) ^log 2 ^»l 

In this limit we expand Eq. (2.4) and Eq. (2.5) around Qo/Q — > and find the leading 



terms 



^.secondary 



^primary 
4 

^secondary 



a, 2 Q / 4 «s ! Q , Qo 
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01 
Q 2 



(2.6) 



Their ratio scales like 0-P nmar yy0- second ary ^ .y^r \ogQ/Q , i.e. the primary emissions 
are logarithmically enhanced. In the limit of a large logarithm (high single emission 
probability) the distribution of final state emissions are increasingly primary, and 
therefore give a Poisson distribution. 

Physically interpreting Eq. (2.6), a second logarithm in the secondary contribution 
would come from the right-most Sudakov of Eq. (2.5). However, it has vanishing 
support for Q — > oo and does not appear in the approximate result. The emitted 
gluon in this case spans a vanishing relative fraction in virtuality space where it may 
emit an additional parton. 

(2) ^log^«l 

Taking this limit of Eq. (2.4) and Eq. (2.5), we find 

r a 2 s 4 _Q 

4(2vr)2 ° g Q 1 ~ V~ s ^ Qo 



a primary(Q 2) Q2) = primary log 4 ^ + Q ( Q 3 log 6 Q 



primary 

= 6 ^d^- SCCODdary W 2 ^o)- (2-7) 

The two contributions become logarithmically equivalent and differ by an 0(1) con- 
stant depending primarily on color factors. In this regime the emission probability 
is small and the final state is selected democratically. The formerly Poisson scaling 
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pattern receives large contributions from subsequent or secondary splittings. Note 
that to justify the logarithmic expansion we still require log 2 Q/Qo > 1 but not large 
enough to spoil the small emission probability. 

In Appendix A we use a toy model of secondary splittings to show that c secondar y gives 
the subsequent splitting parameter n' in an iterated inhomogeneous Poisson distribution. 

2.2 Generating functional for jet fractions 

Following the argument presented in the last section we need a way to derive scaling 
patterns for arbitrarily high parton multiplicities. The generating functional formalism for 
QCD allows us to calculate resummed jet quantities [9, 10]. We construct a generating 
functional in an arbitrary parameter u by demanding that repeated differentiation at u = 
gives exclusive multiplicity distributions. Different moments of the same generating 
functional then produce more inclusive jet observables. This gives us a set of coupled 
integral equations which we can solve in the limit of large and small emission probabilities. 
We will find that the derived jet multiplicity distributions follow a Poisson or staircase 
pattern, respectively. 

For a known series of functions P n we define 



with P, 



n-l 



n=0 



1 d n 
n\ du r ' 



(2i 



u=0 



Note that for generating functionals we always suppress the argument u. In the application 
to gluon emission the explicit factor 1/n! corresponds to the phase space factor for identical 
bosons. The exclusive jet rates P n are defined in Eq. (1.1). In accordance with that 
definition we only count radiated jets, for the generating functional that means we use 
P n _i instead of the usual P n consistently. 



The quark and gluon generating functionals to next-to-leading logarithmic accuracy 



arc 



%{Q 2 ) 



u exp 



u exp 



rQ 2 

/ dtr q (Q 2 ,t)(<z> g (t)-i) 
[ Q dt(r g (Q 2 ,t)(<s> g (t)-i) + T f (t) 



gg) 



(2.9) 



The splitting kernels are defined in Eq. (2.2); gluon splitting to quarks, described by the 
kernel Tj-(i), is suppressed by a power of the logarithm. Assuming Q > Qo, or a high 
emission probability as discussed in Section 2.1, the largest contribution to the t integration 
comes from the region where t ~ Qq and $> qtg (t) ?s ^^(Qq) » u. Both evolution equations 
then read 



u exp 



dtTjiQ^t) (n-l) 



uAj(Q 2 ) 
A j(Q 2 ) u ' 



(2.10) 
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with the Sudakov factor defined in Eq. (2.1). For the jet rates we find a Poisson distribution 

P -A^Mit! nr n> _ I log a, (q 2 ) i ( . 

P n _l - Aj{Q ) r n _ 1 y or R (n+l)/n - • (2.11) 

Again, the jet counting reflects our convention that of n jets only n — 1 are radiated off 
the hard line. The result in Eq. (2.11) reflects the same underlying physics as Eq. (2.6), 
namely a universal logarithmic enhancement of the primary emission over subsequent ones. 
The latter also covers the sub-leading terms to the pure Poisson distribution, and thus 
determines the size of the leading corrections to Eq. (2.11). For a hard quark as well as for 
a hard gluon line this Poisson distribution contains only logarithmically enhanced gluon 
radiation. According to the approximation of Eq. (2.10), any subsequent splitting of the 
radiated gluons is subleading. 

To investigate deviations from this perfect Poisson pattern we first study jet fractions 
P n for up to four emissions in e + e~ collisions. We analytically derive them using the 
generating functional, Eq. (2.8), and show the results in Appendix B. Expanding them to 
0{a s ) and combining this with inclusive unitarity ^ P n = 1 at each fixed order gives the 
double-logarithmically enhanced contribution for n < 5. They show the expected Poisson 
pattern for the abelian terms oc Cf, 

Pn= [A9{ nP ? (j^dtUQ 2 ^)^)) ■ (2-12) 

The additional gluon Sudakov compared to Eq. (2.11) takes into account that the radiated 
gluons do not split in the primary contributions. The non-abelian terms oc Ca do not 
exponentiate with respect to the qq final state. They deviate from Poisson scaling starting 
at two additional jets and provide sensitivity to the triple- gluon vertex [46]. The secondary 
contribution to two-jet emission is color enhanced via Ca/Cf, but ultimately smaller than 
the Poisson term due to the averaging factor over the second splitting function. For the 
leading double logarithm this is a simple suppression factor of 1/6. In Appendix A we 
compare the LL jet rates and a toy model for an iterated Poisson process. 

In addition to purely non-abelian splittings, mixed primary and secondary contribu- 
tions also deviate from the Poisson pattern. This effect we can study in the average jet 
multiplicity in e + e~ collisions. This observable conveniently singles out all emission his- 
tories which are secondary with respect to the core process. This way we generate the 
highest non-abelian terms at each perturbative order, 

(n je ts) = 2 + — % log 2 y cut (2.13) 

7T Z 



a. N 4 C\C^ /«. y C\C F 10 / / as ) - { 12 

1 vr J 322560 V vr / 58060800 V V vr ' 



6 
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The second term gives the Poisson distribution for abelian emissions off the quark line, 
i.e. n as given in Eq. (2.12). All other terms are leading in Ca for a given power of a s , 
i.e. they originate from a single radiated gluon [10]. Mixed terms of order Cp +1 C™ are 
absorbed into purely CpC™ terms through exponentiation. Therefore, the combination 
(ftjets) — 2 — o. s Cf /(2tt) log 2 y cut reflects the non- Poisson nature of the purely secondary 
emission at leading logarithm. 

Going back to the limiting cases introduced in Section 2.1 we also find a recursive 
solution for small emission probabilities when restricting ourselves to pure Yang-Mills the- 
ory. We take Q ~ Qo, but large enough to define logarithmically enhanced terms. The 
generating functional for a theory with only gluons satisfies 



d%(Q 2 



dQ 2 



UasiQ2) (%{Q 2 ) 



6 



' 2vrQ 2 
$ 9 (Q 2 )x 

f 9 (Q 2 ,Q 2 ) ($ 3 (Q 2 )-1) 
after integrating by parts and defining 

f g {Q 2 ,Ql) = Ca 



L) t / dt^- {%{€}-!) 

hi 1 



(2.14) 



C A 
2vrQ 2 



a s (Q 2 
2ttQ 2 



dtlog 4 (< ^ (t) " r 



log 



Q^_ 11 

Ql 6 



(2.15) 



Only including the leading logarithms this is simply the negative splitting function from 
Eq. (2.2), but it obviously differs beyond this approximation. To establish the accuracy of 
Eq. (2.14) including the maximal amount of NLL contributions, we expand in powers of 

^(Q 2 -Qt)/Q 2 , 



r g (Q 2 ,Q 2 ) (<^(Q 2 ) - 1) 

t ( d 



CaOs 
2vrQ 2 

C a ol 



11 



+ 2 



Ca f Q 

^Q 2 U ^Q 2 ^^^- 1 ^ ~ 2ttQ 2 



11a A 
6 ' V" 3vr 

[-4e 2 + C(e 3 )] . 



e + 0(e 2 

(2.16) 



The logarithmic and finite contribution from the first term on the RHS of the differential 
equation are of order e 1 and e° while the formally infinite series from the second term in 
Eq. (2.14) starts at e 2 . Keeping only terms linear in e we obtain the simple form 



d%(Q 2 ) 
dQ 2 



<S> g (Q 2 )r g (Q 2 ,Q 2 )(<t> g (Q 2 )-l) 



(2.17) 



Including the boundary condition $> g (Qo) = u we can solve this 
1 



uA g (Q 2 ) 



with 



A g (Q z ) = exp 



dtt g {t,Q 2 ) 



. (2.18) 
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Figure 2. Jet ratios R( n +i)/n in e + e~ —> jets production at 2 TeV center-of-mass energy. We show 
a Poisson fit with n = 8.7 and a staircase fit to the tail. We use Sherpa [33, 48] with the g — >• qq 
shower splittings switched off. 



Neglecting the effects of the running coupling, A g (Q 2 ) is a Sudakov form factor. Including 
the running coupling, Eq. (2.18) differs from the standard Sudakov in Eq. (2.1) starting at 
higher orders, 



A 9 (Q 2 ) 
A 9 (Q 2 ) 



exp 



o 



12vr 



log 3 % 



(2.19) 



Taking derivatives of the generating functional in Eq. (2.18) at u = we can compute the 
exclusive jet rates 



P, 



n-l 



n-l 



or 



R 



(n+l)/n 



1 - A g (Q 2 



(2.20) 



These constant ratios define a staircase pattern. Comparing Eq. (2.11) and Eq. (2.20) 
we see that in two distinct phase space regimes we find two clear scaling patterns for the 
Yang-Mills or pure gluon case. Both of them can arise in final state gluon radiation, which 
means they should in principle be observable in e + e~ —> jets events. 

The all-order theoretical predictions for Poisson scaling, Eq. (2.11), and staircase scal- 
ing, Eq. (2.20), we can compare to simulated e + e~ — > jets events. To cover both, a large 
scale separation Q S> Qo as well as a democratic scale Q ~ Qo, we use a large center- 
of-mass energy of 2 TeV and a very small lower cutoff y C ut = 5 • 10~ 7 for the Durham 
jet-reconstruction algorithm [47]. In Fig. 2 we show jet ratios Rf n+1 y n for a large range 
of n. Indeed, we observe Poisson as well as staircase scaling. The same behavior is known 
from hadron colliders for example in pp —> 7+jets production [8]: for relatively low n values 
the emission is dominated by large scale differences, inducing a Poisson pattern. For large 
jet multiplicity individual emissions are not affected by a large scale difference, so we see 
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Figure 3. Contribution from secondary emission to the squared matrix elements given in Eq. (2.21) 
as a function of Ar) for = 0, 7r/6, 7r/4, 7r/3 versus the constant primary contribution (horizontal 
line). All curves are normalized to K ~ §^C^/(pt,iPt.2)- 



a staircase tail. While this transition is a solid QCD prediction it has not been studied 
experimentally (yet). 

2.3 Matrix element corrections 

In all of the above discussion we only assume logarithmically induced emission and neglect 
any kind of phase space effects. A simple test case for the relative contributions of primary 
vs subsequent emissions including additional phase space information is two-gluon emission 
from a qq dipole. The squared matrix element for strongly ordered two-gluon emissions 
is [49] 



\M{piiP 2 ) 



32C F 

PT,lPT,2 



C A 



cosh(r/i - 7/ 2 ) 



cosh(r/i - 772) - cos((/>i - 4> 2 ) 



l)+2C F 



(2.21) 



where iji are the gluon rapidities, (pi are the azimuthal angles and pT,i the transverse 
momenta. The term proportional to Cf with its very simple kinematic structure represents 
the primary emissions. The Ca term corresponds to the subsequent emission contribution 
and contains an interesting dependence on the phase space. 

The transverse momentum integration leads to log Q/Qo terms for both contributions, 
which we now assume to not be too large — otherwise we would be logarithmically dom- 
inated and the results from the previous section would apply. Instead, we are interested 
in the 77 and (p dependence. The question is whether such final state kinematics can be 
observed as deviations from our scaling patterns. 

In Fig. 3 we examine the primary and subsequent contribution as a function of A77. 
Testing different A(j) values, we see that the secondary emission dominates when the gluons 
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are close. For A(p = ir/2 the Ca term vanishes for all rapidity separations. Note that for 
y^Af?) 2 + (A0) 2 < R the two emissions are clustered in a single jet of radius R (typically 
R 0.4 — 0.7) and as such do not contribute to the 4- jet rate. This suggests that for 
large angle emissions the nj e ts spectrum will remain Poisson even including matrix element 
information. Or in other words subsequent emissions breaking the Poisson pattern are 
unlikely to be widely separated. 

Although Eq. (2.21) contains particular matrix element information, it is still based 
on the eikonal approximation and does not include kinematical constraints. However, we 
know that subsequent emissions lead to deviations from Poisson scaling in all phase space 
regions. In order to get a handle exclusively on effects from outside the soft or soft-collinear 
regime, we consider the photon multiplicities in the QED process e + e~ — > e + e~ + nj. To 
leading order all photon emissions are primary. 

We consider the cross sections for the processes e + e~ — > e + e~ + fry at y/s = 500 GeV 
using the exact tree-level matrix elements regulated using the kx measure 

2mm(£' 2 , E?) 

Y (1 - cosfltf) > y cu t • (2.22) 

For small values of y cu t we should find a Poisson pattern in the exclusive photon rates, which 
we confirm in Fig. 4. For larger y cut the different multiplicity distributions start deviating 
from the Poisson pattern. The ratios are pushed apart from one another, opposite to what 
we expect from a staircase pattern. The reason is that each emission takes a non-negligible 
amount of the total energy of the event and suppresses the phase space for subsequent 
emissions. Going back to the two main scaling patterns this means that matrix element 
and final-state phase space effects are not responsible for the transition from Poisson to 
staircase scaling. 

3 Hadron colliders 

The analytic form of the fixed energy jet fractions given in Appendix B indicates that a 
final state cascade initiated by a qq pair follows neither a Poisson nor a staircase scaling 
pattern at low multiplicities. However, from Z+jets production we know that essentially all 
jet ratios are constant, with an even stronger suppression of R\/q [50]. This suggests that 
additional effects drive the jet ratios at hadron colliders towards a staircase pattern. One 
possible cause is that incoming partons do not on average carry the same energy fractions x 
for different final state jet multiplicities; in that case we might observe an initial-state phase 
space effect. Second, jets at hadron colliders are typically generated through initial state 
radiation, and we know that the initial-state parton shower behaves somewhat differently 
from final state splittings. 

3.1 Generating functional for incoming hadrons 

The basis of the QCD treatment of hadron collider physics is collinear factorization which 
allows us to employ the generating functional method [9, 10, 14]. Before we can apply 
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Figure 4. Normalized ratios for the photon multiplicity in e + e~ — > e + e~ + 77,7 as a function of the 
resolution parameter y cu t- The solid lines correspond to the perfect Poisson hypothesis. We use 
Sherpa [33]. 

any of this to jet counting we need to clarify our choice of the factorization scale fj,p 
in exclusive raj e t rates, i.e. in the presence of a jet-counting or jet- veto scale py. The 
resummation properties of the DGLAP equation identify the combined renormalization and 
factorization scale with a collinear cutoff below which initial state splittings are unresolved 
and influence only the functional dependence on the partonic energy fraction x. Because we 
are interested in radiated jets with pr > pv we identify the factorization and the jet-veto 
scale, i.e. [if = Pv- Note that this choice furthermore avoids generating additional finite 
though potentially large logarithms in the ratio hf/pv [17]. 

Symbolically, going from final state radiation in e + e~ collisions to deep inelastic scat- 
tering (DIS) with initial state radiation and parton densities we replace the two generating 
functionals, distinguishing time-like from space-like splittings, 

<*> q (Q 2 ,Pv) x $g(Q 2 ,Pv) -> 3,/ ? (Q 2 ,j£)x Z q/q {x,Q 2 ,Pv)- (3-1) 

As in Sec. 2.2 we omit the argument u in all generating functionals. In the original DIS 
context all scales are defined in terms of the e + e~ Durham algorithm [14], most notably 
the hard scale Q and py = (j,p as well as the softer resolution scale Qo < (ip- We identify 
all three relevant scales Qo = \ip = py. For the DIS analysis this corresponds to not 
further resolving the original macro-jets which define the separation of resolved jets and 
beam jets [14]. Again, this choice omits potentially large finite scale logarithms in our 
perturbative treatment. 

We also introduce an explicit x dependence in the generating functional for incoming 
partons as it is clear that PDF effects alter the possibility to radiate jets. Each emission 
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takes away an energy fraction 1 — z of the emitter; the x value has to change correspond- 
ingly and splitting between different partons needs to be taken into account. From the 
factorization theorem we know that PDFs and partonic cross-sections also factorize at the 
generating functional level, 

l 

Z a (x,Q 2 ,p 2 v ) = J2 J~ f b ( X z ,P 2 v)z b a (z,Q 2 ,P 2 v)- (3-2) 

b x 

The parton densities we consistently evaluate at the scale py. This way, logarithmically 
enhanced parton splittings above py are described by the partonic generating functional 
Z b a . For the generating functional in DIS we start with a time- like generating functional 
for a single (anti-)quark and weight it with the proper electromagnetic coupling [14] 

$dis = $> a {Q\p 2 v ) Z a (x,Q 2 ,p 2 v ). (3.3) 

a 

The partonic cross-sections and jet evolution are the same for the quark and the anti-quark, 
but the PDFs are different. In DIS the final state kinematics fix x. Additional jets radiated 
off the incoming parton imply that in this case we probe higher x values as given by the 
convolution in Eq. (3.2). 

The task is to find the evolution equations for Z\. To leading logarithm (LL) this turns 
out to be relatively simple. In the soft and collinear limit [9, 10] the eikonal approximation 
implies z ~ 1. Furthermore, the g — > qq splitting is logarithmically suppressed compared 
to the other splittings, so we can neglect it. Under these two LL assumptions the evolution 
equation and the corresponding generating function in Eq. (3.2) read [14] 

Z b a (z, Q\py) = 5(1 - Z) 5 b a <S>a(Q 2 ,P V ) 

Z a (x,Q 2 , P 2 y) = fa(x,p 2 y) MQ 2 ,P V ) • (3-4) 

The PDF effects and the jet generating function factorize in x, so we can treat them 
independently. In general, <3? a is a two-scale generating functional [14, 51]. Because we 
identify Qo = py the second scale is suppressed and its evolution equation is almost the 
same as in the e + e~ case. To leading logarithm we find 



^a(Q 2 ,Py) =exp 



I dtT a (Q 2 ,t) ($ g (t,p 2 v )-l) . (3.5) 

Compared to Eq. (2.9) the factor u in front of the exponential is missing. The reason 
is that we cannot resolve a jet if there is not at least one space like splitting. The hard 
parton cannot produce a final state jet, so we always find the normalization condition 
& a (Py,Py) = 1- The further evolution of emitted partons we describe with the time-like 
functional of Eq. (2.9). 

Moving on to Drell-Yan production with two incoming partons we need to replace the 
generating functionals, symbolically written, to 

Z q /q(Xa,Q 2 ,Py) X Zg /q (x b , Q 2 , Py) . (3.6) 
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Thus, we replace the remaining time-like generating functional with a space-like generating 
functional to describe two incoming partons. A major complication is that the final state 
phase space does not fix x a ^ anymore. Instead, we have to integrate over their allowed 
ranges and find the generating functional for the Drell-Yan process, 




~ £/a(* (0) ,^)$a(Q 2 ,P^) MX {0 \P 2 V )MQ 2 ,PV) • (3-7) 
a, b 

From this generating functional we can derive the individual n-jet rates. For the second 
line of Eq. (3.7) we use the leading logarithmic approximation as in the DIS case. To 
arrive at the third line we replace the variable x a ^ values by a typical partonic energy scale 
x( '. For typical hadron colliders processes we assume this value to be close to threshold 
and equal for the two incoming partons. The argument u which generates the different 
n-jet rates is carried only by the generating functionals & a ,b(Q 2 \P\/)- Starting with two 
generating functionals for the two initial state particles, hard jet radiation with px > Pv 
indeed factorizes from a PDF factor. 

One apparent contradiction related to the PDF kinematics we need to resolve. On 
the one hand, in Eq. (3.4) the eikonal approximation allows us to set z « 1, which means 
that the entire energy dependence is encoded in the PDF factor. On the other hand, each 
resolved jet requires a finite px > Pv- Hence, the integration range for x a ^ is determined 
by the partonic n-jet process and x^ implicitly depends on u. This implicit dependence 
we have to account for by hand. In particular for parton density regimes which increase 
towards small x the majority of multi-jet events at the LHC are produced at threshold. 
The threshold value for any of the n-jet production rates we denote as x^ n \ leading us to 
the modified factorized form 

^Drell-Yan = £ fa (x {n) , P V )<5>a (Q 2 , Py) fb(x {n \ p 2 V )MQ 2 , Pv) ■ (3-8) 
0,6 

We emphasize that n is determined a posteriori upon differentiation with respect to u, so 
is presented for illustrative purposes only. Eq. (3.8) means that to leading logarithm the 
jet radiation pattern in the Drell-Yan case is the same as in e + e~ — > jets processes, modulo 
explicit PDF factors estimated using an n-dependent threshold kinematics. A similar 
approach can account for energy momentum conservation in soft-gluon resummation [52]. 
This way we leave the LL evolution of jets untouched and instead shift the x value in the 
PDFs to account for additional jets. All our findings from Sec. 2 we can immediately apply, 
once we understand the PDF correction factor in the next section. 

3.2 Parton density suppression 

In Sec. 3.1 we have learned that to leading logarithmic accuracy the effects of the parton 
densities and jet emission factorize. For large jet multiplicities this explains the observed 
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staircase scaling at hadron colliders [7, 41]. Parton densities contribute to this effect in 
particular at low multiplicities. When increasing the jet multiplicity the typical partonic 
energy fractions x probed by the partonic process increase as well. The relative increase 
in x is largest for low jet multiplicities. 

In terms of the assumed threshold kinematics adding a jet with finite transverse mo- 
mentum implies x( n+1 ) > x^ n \ To compute the relative cost of producing an additional 
jet we estimate the ratio of PDF values evaluated at x^ and x^ n+1 ^ as a function of the 
number of extra jets n. In effect this is the discretized second derivative with respect to 
x. For hadron collider processes involving two parton densities f(x, Q) we define the PDF 
correction factor to the ratio of successive jet ratios R( n +i)/n/ R(n+2)/(n+i) 



/(a (w+1) ,Q) 



/(x (n+2) ,Q) 

/(x(»+D,Q) 



(3.9) 



The square in the definition of B n reflects the two PDFs in hadron collisions. If for example 
the partonic ratio of two successive jet ratios is R( n +i)/n/ R(n+2)/(n+\) ~ c then the proper 
hadronic ratio becomes B n c. We fix Q for simplicity, but this only mildly affects our results. 

The main effects are, first, that B n < 1 in most cases. This way PDF effects suppress 
the lower multiplicity ratios R( n +i)/n- For large jet multiplicities the relative impact of 
yet another jet becomes small, B n — > 1. The hadronic initial-state effect on the jet scaling 
disappears and we are back to the staircase pattern. Second, the PDF effect is largest 
for the steep gluon densities, as compared to the flatter valence quarks. Finally, allowing 
for variable Q the PDF values f(x, Q) increase (decrease) with Q for low (high) x, with 
a cross-over point around x ~ 0.1. For small x values the initial state evolution then 
suppresses jet ratios at high multiplicity or large Q 2 . 

What we are most interested in are PDF effects for the Drell-Yan process at lower 
multiplicities. We consider the threshold values x^ n \ for example for producing an on-shell 
Z-boson and one additional jet, 



<m% + 2 (pt\ Pt + m 2 z +p T ) 
x« = -V 1 . (3.10) 

2 -Kbeam 

where £ beam = 3500 GeV for the LHC in 2011. Comparing ajM with a;(°) « m z /(2£ beam ) 
shows a sizeable shift. For the two-jet threshold x&> two limiting cases are either the addi- 
tional jet adding merely PT/(2-Ebeam) to x^- 1 ' or two approximately collinear jets recoiling 
against a hard Z. The variation between these two cases estimates the uncertainty on our 
method which can be generalized straightforwardly to the n-jet final state. 

In the left panel of Fig. 5 we display B n for the estimated Drell-Yan kinematics, 
assuming each jet has transverse momentum Pt = Pv = 30 GeV. The effect on the first jet 
ratios is large, but quickly diminishes towards higher n. We also see that if we require a 
leading jet with large transverse momentum, pl^ ad > 100 GeV, we move to sufficiently high 
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Figure 5. Left panel: estimated PDF suppression for inclusive (solid) and jet-associated (dashed, 
pi^ad > GeV) Drell-Yan kinematics. We assume an initial state with d-quarks only. Right panel: 
same for Higgs production in gluon fusion with mn = 125 GeV. The uncertainty encompasses two 
representative kinematical limits of the multi-jet final state, described in the text. 



x such that additional jet ratios are unaffected by the PDF effect. It is reassuring to see 
that if we combine the PDF suppression of R\/q (0.46 - 0.65) with the Ca enhancement of 
i? 2 /i (1.36) and assume an original Poisson scaling we find R1/0/R2/1 = (0.67 — 0.95), in 
nice agreement with ATLAS data [50]. This beautifully illustrates that staircase scaling at 
large multiplicities can be derived from first principles QCD while for small multiplicities 
it is something like a sweet spot. 

As an additional check, we present the PDF suppression in gluon-fusion Higgs produc- 
tion in the right panel of Fig. 5. We assume m# = 125 GeV, ignore flavor changes and 
consider jets with = pv = 30 GeV. The gluon PDF drops more rapidly for increasing 
x, inducing a large PDF suppression. On the other hand, the increasing energy of the core 
process as compared to the Drell-Yan process slightly decreases the effect. The combination 
of the two gives remarkably similar results to the Drell-Yan process. 

3.3 Initial-state parton shower 

As indicated above, jet radiation at hadron colliders is generated mostly through initial 
state radiation, which means that our final-state analysis of Sec. 2 should be modified. We 
need to compute the spectrum of jets arising as primary emission in the backward evolution 
and acting as seeds for subsequent final-state radiation. For simplicity we just consider the 
backward evolution along an initial-state quark line with a single type of branching, namely 
gluon emission. The evolution proceeds through the space-like Sudakov form factor 

* 2 dt f dz a s n ( A f q (x/z,t) 



n(t 1 ,t 2 ;x) = ex P p / y j 7 ;W^h (3.11) 
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Figure 6. Normalized ratios of gluon production rates originating from the backward evolution 
of d-quarks according to Eq. (3.11). We assume different initial values for x^°\ Each splitting is 
restricted to z values generating a minimal increase of Ax = 0.02. The dashed line indicates the 
expectation for a perfect Poisson distribution of the underlying jet rates. 

cf. Ref. [53], with the appropriate splitting kernels for gluon emission Pq — >qg and ignoring 
potential initial-state flavor changes here. The evolution of each initial state parton starts 
with momentum fraction Xj, determined by the hard process, and virtuality t = X1X2S, 
it terminates at the hadronic scales x ~ 1 and t = Qq associated with the transition to 
non-perturbative physics. 

We know from Sec. 2 that a single time-like Sudakov form factor produces perfect 
Poisson scaling. Indeed, whenever we have a non-emission probability represented by e _r 
as in Eq. (2.3), where T does not change as a result of a splitting, the process is guaranteed 
to produce Poisson scaling. For backwards evolution the situation, however, is different. 
Once a splitting is generated, e.g. using a veto algorithm [54], we need to re-compute 
x — ¥ x/z because each emitted parton increases the combined x value. The evolution then 
proceeds with this different effective splitting kernel. In other words, the PDF dependence 
in Eq. (3.11) explicitly correlates parton emissions, breaking a key ingredient to the Poisson 
derivation. 

To quantify this effect we numerically evaluate the gluon emission spectrum generated 
by the Sudakov form factor given in Eq. (3.11) using a veto algorithm. We neglect any 
recoil effects and therefore expect a somewhat smaller suppression in the first bin. In Fig. 6 
we display the normalized ratios of the gluon production rates off initial state <i-quarks. 
We assume different starting values for x*- ) thereby keeping the evolution distance and the 
minimum step size of Ax = 0.02 fixed. This roughly corresponds to the emission of a hard 
additional jet with pt ~ 70 GeV which we are interested in for the LHC. The starting scale 
we vary between x^ ^ = 0.03 and x^ ' = 0.2 corresponding to a hard process of 105 GeV 
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and 700 GeV respectively. The jet ratios strongly deviate from the Poisson pattern for low 
x(°\ while the effect quickly diminishes for large x^°\ This is the same pattern we find for 
the PDF effect in Sec. 3.2. 

Finally, to see how our two approaches to parton density effects are related we study 
the PDF part of the weight attached to a single resolvable emission coming from a collision 
with momentum fraction xq. Using the backward evolution Eq. (3.11) the exclusive one-jet 
rate can be represented as 

dt f dz . f(xW/z,t) 



*i ~h(x {0) ,Q 2 )f2(x^,Q 2 )ar OmC fjj d -^P{zY 



/(*(«), t) 



f f~(0) n 2 \ f r~(0) n 2 \ ^partonic f{x^\Q 2 ) P(z^) 

~ h(? y > ,Q ) h(? y > ,Q ) ° f( x (°),Q 2 ) ' ( ' 

In the second line we limit ourselves to the leading logarithmic approximation (i.e. ignoring 
the t dependence of the PDFs) and fix the resolvable momentum fraction to its threshold 
value z^ 1 ' = x^/x^\ The effect of the PDF weight in the ratio cti/oo then turns into a 
suppression factor /(x«,Q 2 )//(x(°),Q 2 ). Each emission just pushes up the overall PDF 
suppression and we are effectively led to the estimate on the shape of the ratios provided 
by Eq. (3.9). 

3.4 BFKL evolution 

So far, we only consider the DGLAP evolution which relies on collinear factorization and 
resums collinear logarithms. Parton evolution can also be represented by BFKL [55] or 
CCFM [56] dynamics which rely on an entirely different form of the factorized matrix 
element. In this approach a simple expression for the n-jet generating function at leading 
logarithmic order in logl/x and log Q/pv reads [57] 



®{Q ,£V)bfkl = exp log — 

TTW pv 



1 + (1 - u) log — 

TTW pv 



u/(l-u) 

. (3.13) 



In this expression w is the Mellin conjugate variable of x which for the physical jet rates 
requires convolution with the structure functions and transformation back to x space. In 
analogy to Eq. (2.9), functional derivatives of <J?bfkl evaluated at u = return exclusive 
jet rates. 

As in the DGLAP case of Sec. 2.1 we can compute the scaling patterns in the limit of 
small and large emission probability. 

(1 log — > 1 

TTW pv 

Following Eq. (3.13) this large logarithm describes the limit of large emission proba- 
bilities. Again, we find that the nj e ts distribution shows a Poisson scaling 

u/(l-u) 



1 d n 



1 + (1 - u) log — 

TTW pv 



— log 1 + log — 

n! \ ttw pv 



u=0 

(3.14) 
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Note that this result is obtained by taking the limit after the differentiation. Although 
Eq. (3.14) is formally true in the limit \jw log Q/pv oo the enhancement of the 
Poisson term is gradual, i.e. log n (l/u; log Q/pv)- At every multiplicity there appear 
terms of the order log^ 1 (1 /w log Q/pv) with possibly large coefficients. Therefore, 
we expect this Poisson distribution to only emerge at very high energies and currently 
experimentally inaccessible x values. 

(2) 2 ^log^«l 

Expanding to leading order in the emission probability the rates are 

P. - (^logQ)". (3.5, 

\ 7TW pv ) 

This is a staircase distribution in the jet ratios. This result can also be seen in the 
corresponding fixed order computation [58]. Apparently, this staircase distribution is 
of an entirely different origin from the gluonic cascade in DGLAP, Eq. (2.20). In the 
BFKL evolution there is no notion of subsequent emissions, all gluons are emitted 
directly from the factorized i-channel. A subsequent emission does not contain the 
logl/x enhancement. For the physical jet rates, we transform to x space via a 
convolution and the rates are of course different; only the scaling, i.e. the shape of 
the ratio distribution, does not change. 



Unlike the patterns we see in the DGLAP approach, both limits pointed to above 
are only realized in corners of phase space which are not particularly relevant to the hard 
processes at the LHC. The actual nj e ts distributions in BFKL- type regions will not be 
described in terms of a pure scaling pattern. In fact, it was pointed out in Ref. [59] that 
the njets distribution from simulated BFKL emissions are described better by a shifted 
Gaussian than a Poisson distribution. From a scaling perspective it is difficult to imagine 
how this might be induced by first principles QCD. 



4 LHC predictions 

ATLAS and CMS analyses based on 2011 and 2012 data with an excellent understanding 
of the detectors allows for detailed investigations of high-jet-multiplicity final states. The 
scaling patterns derived in this work should hence be experimentally testable before the 
2013 shutdown. 

From the discussion in Sec. 3 we know that n-jet scaling arguments can be transferred 
from e + e~ to hadron colliders without major modifications. This means we can search 
for Poisson scaling and staircase scaling in standard candle processes at the LHC. For 
reasonably large jet multiplicities both of these patterns are derived from first principles 
QCD. Parton density effects modify the first few bins in the exclusive jet ratios R( n +i)/n- 

In the definition of the observables we need to carefully distinguish between inclusive 
and exclusive jet rates, i.e. counting exactly n jets or at least n jets in the n-jet bin 
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radiated off a given hard processes. Staircase scaling is unique in the sense that it appears 
in exclusive and inclusive observables with the same constant ratio R( n +i)/n [!]• Poisson 
scaling at e + e~ colliders or hadron colliders is limited to exclusive jet rates. Its translation 
into inclusive observables is tedious [1] . 

Following Sec. 3 we expect that QCD features together with the PDF effects produce 
a convincing staircase pattern for cross sections with a democratic jet selection and in the 
absence of major kinematic cuts [1, 8]. The ratios R( n +i)/n should be constant, Eq. (1.2). 
Towards large jet multiplicities this constant value R is largely independent of the under- 
lying hard process. The first nj e t s bin is most sensitive to the parton densities, as shown in 
Sec. 3.2. The PDF effect will generally lead to a suppression of the otherwise large value 
of Rx/o, extending the staircase pattern to low multiplicities where the QCD derivation 
fails. However, the PDF effect is clearly not independent of the hard core process which 
can involve incoming gluons as well as incoming quarks. 

In contrast, for a jet selection with a large ratio between the core-process scale and the 
jet acceptance cut {pr,j > Pv) we expect several bins of the nj e ts distribution to follow a 
Poisson distribution. Only at large jet multiplicities the ratios will overshoot the Poisson 
pattern n/{n + 1), turning into a staircase tail. This is the same behavior we derive for 
e + e~ — > jets production in Sec. 2. 

Our considerations about the nj e ts distributions associated with a general hard process 
needs to be validated for a variety of channels, potentially in dedicated analyses. As a 
starting point, we propose two simple extensions of already existing analyses. They serve 
as stringent tests in particular of the Poisson hypothesis for the low-multiplicity bins and 
the transition to a staircase-like behavior for higher nj e t s bins. We start with the Drell-Yan 
process, but asking for a leading jet with significant transverse momentum p^ ad > 5py 
while all additional jets are selected democratically with px > Pv- Our second example is 
the njets distribution inside a rapidity interval defined by a dijet core system. 

4.1 Z+jets 

For the VF/Z+jets channels recent LHC measurements of the nj e t s distribution support 
a staircase hypothesis [41, 50]. The key ingredient is the democratic selection of jets, 
e.g. pr,j > 30 GeV for all jets in Z — > u + u~ and Z — > e + e~ events [50]. We propose 
a minimal modification of this analysis by increasing the pt cut on the leading jet to 
pkad y 150 GeV. The hard core process is now defined as the Z boson combined with the 
leading jet, the first of nj e t s radiated QCD jets is the second jet in the event. According to 
Sec. 3.2 such a selection makes the PDF suppression effect rather marginal and we expect 
to see a Poisson distribution in the low- multiplicity bins. 

To simulate this process we use Sherpa vl.4.0 [33, 48]. We employ Sherpa's tree-level 
merging algorithm based on truncated showers [61], including real-emission matrix elements 
for up to five final state partons, the merging scale we set to 25 GeV. In Fig. 7 we present 
the exclusive nj ets cross section ratios, keeping in mind that the counting of the radiated 
jets does not include the hard leading jet which together with the Z boson constitutes the 
core process. As expected, we observe a clear Poisson scaling for nj ets < 3. This can be 
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Figure 7. Left panel: exclusive jet ratios for Z + jets production at vS = 7 TeV. We require a 
leading jet with p^ ad > 150 GeV. All other jets have Pt > Pv = 30 GeV. The line shows a Poisson 
shape with n extracted from the first bin. Right panel: cross sections as a function of the event-wise 
ratios of the jet transverse momenta in inclusive Z + 3 jet production. While the leading jet has to 
pass p^ ad > 150 GeV additional jets are selected uniformly with just pt > Pv = 30 GeV. 



seen when comparing the individual bins to the Poisson shape R( n+ iy n = n/{n + 1) where 
n is fixed by the first bin Ri/q. 

However, already the three-jet rate comes out higher than the Poisson extrapolation. 
At this point the non-abelian nature of QCD radiation takes over, giving us a staircase tail 
with R ~ 0.36. This behavior is strongly reminiscent of Fig. 2 for e + e~ — > jets. An updated 
analysis of the 2011 data sample and certainly the 2012 data set should allow for a test of 
up to five additional jets, probing QCD predictions for the complete rej e t s distribution. 

In addition to its defining power of very general features of the nj e ts distribution the 
jet selection cuts are reflected in the transverse momentum distributions of the jets. In the 
right panel of Fig. 7 we depict the cross sections as a function of the ratios of the py-ordered 
jet transverse momenta defined event-by-event. For this distribution we for once deviate 
from our usual exclusive jet counting and consider events with at least three jets in the 
final state. We consider the pt ratios first-over- leading (black), second-over-leading (green) 
and second-over- first jet (red). With the leading jet px > 150 GeV and the additional jets 
selected uniformly with pt > 30 GeV we expect the first and second jet to peak around 
the selection cut. This is confirmed by the simulated results that exhibit strong peaks 
for the corresponding ratios around p^ st /pJ^ ad ~ 0.2 and p^? d /plf, ad ~ 0.2. However, it is 
interesting to note that the ratio p^ id /p^f t does not peak around 1. Rather QCD favors 
the first radiated jet to be significantly harder than the second. The obtained distribution 
in fact turns out to be more or less flat between 0.25 and 1. It is certainly interesting to 
study these observables in addition to the raj e ts distribution, as they contain complementary 
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Figure 8. Left: Predictions vs ATLAS data [65] for the di-jet gap fraction (py = 20 GeV) as a 
function of Ay and in slices of Pt- Right: average number of gap jets as a function of Pt in slices 
of Ay. The tagging jets are defined as most forward and most backward. All jets are reconstructed 
using the anti-fc t algorithm with R = 0.6. 

information on the underlying QCD dynamics. 
4.2 QCD gap jets 

An interesting set of observables from the perspective of multi-jet final states are gap frac- 
tions or gap jets. In that case we require a specific kinematic structure of hard and widely 
separated jets and count the QCD jets in between. The core process is the production of 
two widely separated hard jets. 

A recent ATLAS study [65] identifies two forward jets, so-called tagging jets, either as 
the highest px (selection A) or the most forward and backward in rapidity (selection B). 
The core di-jet system is defined in terms of px = (pr,i + Pr,2)/2 and Ay = \yi — 7/2 1 - The 
gap fraction Pq = <7o/crt t ^ s given by ai l events with no additional jet in between the two 
tagging jets and px > Pv = 20 GeV. ATLAS measures it as a function of Ay. In addition, 
they provide us with the average multiplicity of vetoed jets or gap jets. Following Sec. 2 
a hard cut on px will typically enforce a larger logarithm than simply log(px/pv)- This is 
because the two tagging jets tend to be asymmetric in px- As a consequence we expect a 
Poisson scaling for the gap jets. 

Perturbatively, the gap fractions are dominated by the single emission probability 
while the average number of gap jets probes multiple emissions. We again propose a 
simple extension of the existing analysis. In addition to the average number of gap jets 
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Figure 9. Jet ratios for exclusive gap jets in slices of pr for the tagging jets, denned as the most 
forward and backward in the event. The rapidity separation of the gap-defining jets is 1 < Ay < 2 
(left) and 2 < Ay < 3 (right). 



the corresponding exclusive nj e ts distribution should be studied. Reproducing the average 
number of gap jets cannot validate a full nj e t s distribution, where the latter directly tests 
QCD scaling patterns. 

Again, we use Sherpa with tree-level merging, including hadronization and underlying 
event. We account for up to four hard final state partons in the veto region, the merging 
scale set to Q cu t = 20 GeV. This way, the emission of wide-angle soft gluons should be 
correctly modeled including the full color structure. The DGLAP based parton shower 
resums logarithms of large px ratios only, for sufficiently large rapidity gaps one may 
expect the agreement to suffer [66]. In Fig. 8 we compare Sherpa predictions with ATLAS 
data [65] for both the gap fraction and the average number of jets in selection B. For both 
observables the agreement with data is excellent. Even for large Ay > 4 the simulation 
works well within the statistical limitations, so we do not see the merging of tree-level 
matrix elements with the parton shower breaking down. A detailed analysis of the same 
data using the MC@NLO implementation in Sherpa and assessing both perturbative and 
non-perturbative theoretical uncertainties is presented in Ref. [67]. 

With the confidence gained by correctly modeling the gap fraction and the average 
number of gap jets, we proceed to predict the nj e ts distribution for gap jets. In the forward- 
backward selection this includes all resolved jets with px > py. In Fig. 9 we present the 
gap jet multiplicity ratios for two different rapidity separations, 1 < Ay < 2 (left) and 
2 < Ay < 3 (right), and in slices of Pt- Enforcing a large ratio between the px of the 
tagging jets versus the jet-counting scale py = 20 GeV induces a Poisson distribution for 
all the considered px selections. Comparing the two panels of Fig. 9 we see that the shape, 
though not the normalization, of the nj e t s distributions is largely independent of Ay. For 
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large Pt the highest jet multiplicities suffer from non-negligible statistical fluctuations, in 
particular for the 4/3— jet bin. 

For kinematical selections which induce a Poisson scaling we can approximately derive 
the average number of gap jets from Ri/q. This number determines n and allows for a 
straightforward calculation of the veto survival probability or gap fraction, Pq = exp(— n). 
Comparing this prediction to the explicit n measurement of in Fig. 8 we find both values 
consistent. Such a comparison will provide insight into jet-veto survival probabilities, not 
only for the analysis presented here but also in the context of Higgs analyses [1, 68]. 

5 Conclusions 

While it has been known for a long time that the exclusive nj e ts distribution at hadron 
colliders follows simple scaling features [7], a proper understanding starting from first 
principles QCD has been missing. The two underlying scaling patterns for jet emission is 
a Poisson shape or staircase scaling. 

1. Poisson scaling is well known from multi-jet final states in e + e~ production. It can 
be easily derived in the abelian limit of QCD, i.e. gluon radiation off hard quark legs, 
in complete analogy to soft photon emission in QED. Using the parton shower picture 
as well as generating functionals we have found that a Poisson distribution is tied to 
a large hierarchy of scales and a logarithmically enhanced radiation matrix element. 
Such different scales can be induced by kinematic cuts, for example in Z+jets [1] or 
7+jets [8] production. Many effects, like non-abelian gluon splittings, phase space 
limitations, or sub-leading logarithms modify the pure Poisson shape of the exclusive 
njets distribution. 

2. Staircase scaling [7] is defined as constant exclusive jet ratios R( n +i)/n = &n+i/&n = 
R. It has been observed at hadron colliders since UA1 and generally appears for 
a democratic jet selection, avoiding large scale separations. This includes Higgs 
production in gluon fusion and in weak boson fusion [1] . For large jet multiplicities 
we have found that staircase scaling can be predicted as a non-abelian QCD effect. At 
e + e~ colliders as well as at hadron colliders it describes the exclusive nj e t s distribution 
for large jet multiplicities. At lower jet multiplicities parton density effects violate our 
picture of jet radiation off an unchanged hard process. Numerically, this correction 
restores the staircase pattern, often down to the first emitted jet. 

After this derivation and study of the two generic patterns which are based only on 
QCD and patterns of the parton densities this prediction should be tested experimentally. 
We present two minor modifications of current LHC analyses which would be well suited 
to probe the universal nature of our findings. First, the classical Z+jets channel should 
be used to study the transition from Poisson to staircase scaling using a simple kinematic 
cut on the leading jet. A similar proposal exists for 7+jets production [8]. In addition, 
QCD gap jets already probe these scaling features through the measured average number 
of gap jets. Extending the main observable from the mean number of g£tp jets to the Wjets 
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distributions would allow us to confirm what our first estimate suggests — that gap jets 
can be described precisely using the transition from staircase to Poisson scaling. 

While from a perturbative QCD point of view the exclusive number of jets nj e ts is 
not a particularly attractive observable, it clearly has huge benefits in LHC analyses. In 
the future, jet vetoes in Higgs searches [1] are going to become more and more relevant, 
separating the different underlying production processes. In this paper we have shown that 
it is possible to predict the main features of the nj e ts distribution from QCD. Testing the 
universality of these features would give us a new and improved tool for a huge number of 
LHC analyses to come. 
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A Iterated Poisson process 



The process generated by summing an arbitrary number of independent Poisson processes 
P = P{n\) + P{n2) + • • • P(n n ) is itself a Poisson process with the expectation value 

n = ni + rti H n n . In order to deviate from a Poisson shape, we allow for the occurrence 

of subsequent splittings. A QCD example would be gluon radiation off a quark line with a 
subsequent splitting of the gluon into two gluons or two quarks. 

Suppose we start with a single mother Poisson process described by the parameter n. 
Every emission by the mother process generates a daughter process which itself evolves 
through the same scale as an independent Poisson process with expectation n' . In general 
we will have n 7^ n' . For the mother (daughter) process with Poisson parameter n (n'), the 
two parameter distribution can be approximated by 



P(n;n,n') = e'^'T ( ^ ^ 



n ll n n ~ l 



i=0 



HG [1-n, 2, -n/n'] | , (A.l) 



nl 

where HG is the confluent hypergeometric function with integral representation 

HGk b, z] = — — / e-T-Hl + t) b - a ~ l dt . (A.2) 
r («) Jo 

Let us note some observations regarding the physics content of this model: 

1. The exponential Poisson model is reproduced in the limit n' — )■ 0, where the additional 
splitting probability vanishes. 

2. The zero emission probability is unchanged, while the one emission rate includes the 
additional non-splitting probability of the daughter process. It contains an additional 
suppression e~ n ' which universally reduces the first n-jet ratio for the iterated Poisson 
process. 

3. In the limit of a large additional splitting probability n' 3> n the theory displays 
staircase scaling. The highest order term in Eq. (A.l) carries a compensating factor 
of n! in the numerator. 

4. The convergence of Eq. (A.l) is guaranteed for all n and n' . However, inclusive 
unitarity of the cross-section is not possible; the sum over n in Eq. (A.l) is smaller 
than one for all n' 7^ 0. We expect then that our toy model undershoots the full 
calculation, as confirmed by the e + e~ jet fractions (see Fig. 11). 

5. We emphasize that the purely probabilistic model given by Eq. (A.l) does not capture 
all effects of leading- log QCD. In order to promote our toy model to full double leading 
logarithmically accurate QCD as far as the jet rates are concerned, we need to make 
two modifications. First of all we need to include the full Poisson history by hand 
(i.e. no recursive formula is achievable), and second, we require n' to change as a 
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Figure 10. Ratios from an iterated (normal) Poisson process for n = 0.3 and n' = 3 (n' = 0) given 
in Eq. (A.l). In the iterated case we see a general flattening of the distribution. 

function of the "generation" of subsequent emissions. More precisely the fully non- 
abelian ~ CfC\ three-gluon rate is correctly described if n" = (6/5)n'. This value 
for n" then correctly describes the CpC\ coefficient to the four-gluon rate. In this 
manner we can essentially boot-strap the leading coefficients. 

The important consequence of Eq. (A.l) for realistic values of n' > is that the jet 
ratios Rt n +i)/ n are flatter than the naive Poisson expectation. In Fig. 10, we display the jet 
ratios generated by Eq. (A.l) for n = n' . By making analytic estimates for the subsequent 
splitting parameter we can compare our toy model with jet rates predictions from QCD. 

We assume the relation between the two expectation values n and n' to scale with the 
relative splitting probability for the simplest non-abelian emission. The 2-jet fraction at 
0(a 2 ) in double logarithmic approximation then gives n/n' = 12Cf/Ca- For n < 1, the 
specific value of the two Poisson parameters do not affect the shape as long as their ratio 
is held fixed. Therefore, we compare the shapes of the iterated Poisson model with the jet 
fractions by normalizing the first ratio. This also allows a simple comparison to a purely 
Poisson process with ratios l/(n + 1). The results are shown in Fig. 11, where we find 
that the P2/P1 ratios are in exact agreement. Using Eq. (A.l) the ratio is simply n/2 + n' , 
which our choice of n' reproduces precisely. The higher bins slightly undershoot the jet 
fractions in the iterated Poisson model. The reason stated previously is that the n' should 
change as a function of the "generations" and the certain contribution need to be added by 
hand. The first evidence of this is contained at 0(a 3 ) in P3. The jet fraction coefficients 
for the fin' and n' 2 terms are 1/12 and 1/90 respectively, while the iterated Poisson model 
gives 1/12 and 1/144. For higher terms still the discrepancy persists. The conclusion from 
this plot is that the purely final state gluon cascade generates a staircase scaling pattern at 
higher multiplicities but seemingly fails to explain the suppression on the lower multiplicity 
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Figure 11. Ratios of e + e —> qq + ng jet fractions (black) as computed in Appendix B compared 
to the iterated Poisson model with n' = l^Cp/C^n (red), and Poisson scaling (green). 



bins. For this we need to account for PDFs. 



B Jet fractions 

We briefly quote e + e~ — )• qq + ng jet fractions at C(«gL 10 ) in double logarithmic approx- 
imation and for the Durham algorithm. These are used to compare the relative size of 
subsequent terms in a purely final state cascade picture where the single emission prob- 
ability is still smaller than unity. Results can be obtained from Refs. [47, 69-71]. The 
abelian contributions are shown here to all orders such that in the limit Ca — > we recover 
Poisson scaling. Defining a = a s /ir, L = log7/ cu t and Pi = ai/a to t we have: 
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